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1 Introduction 

One of the oldest questions in strong interaction physics is the density of hadrons in the 
spectrum as a function of mass for large mass. It was conjectured long ago by Hagedorn 
that this density (when suitably averaged) grew exponentially with the mass [1, 2]. A 
useful way to parameterize this is via its integral, N(m), the number of hadrons with mass 
less than m. One way to state Hagedorn's conjecture is that at asymptotically large m, 



where Th , the so-called Hagedorn temperature, is a parameter controlling the exponential 
growth. Note that in a simple model, where hadrons are treated as a noninteracting 
free gas, the Hagedorn temperature represents an upper bound on the temperature of a 
hadronic phase of matter as the energy density diverges for T > Th- The power-law 
prefactor plays an important role in attempts to fit the Hagedorn spectrum from data [3] 
and also determines the thermodynamic behavior of strongly interacting matter as Th is 
approached from below in the simple noninteracting hadron gas model [4]. A more useful 
way to state Hagedorn's conjecture for the purposes of this paper is that for asymptotically 
large masses there exists positive value of T such that 



Th is the maximum value of T for which eq. (1.2) holds. 

As shown in figure 1, the extracted masses of hadronic resonances [5], N{m) does, 
indeed, grow very rapidly up to the point where it becomes difficult to extract resonance 
parameters from experimental data (around 2 GeV). This behavior appears to be consistent 





N{m) > e m ' T ; 
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Figure 1. N(m) for nonstrange mesons using mesons masses extracted from various hadronic 
processes extracted reported by the Particle Data Group [5]. The fit is of the form N(m) = 
am b e m / TH and yields a Hagcdorn temperature of 426 MeV. The fits were done for mesons with 
masses up to 2300 MeV. Our estimate of the Hagedorn temperature is consistent with the results of 
ref. [3]. Note that it is almost 2.5 times larger than the critical temperature Tq obtained by lattice 
gauge calculations [6]. 



with the notion that QCD does have a Hagedorn spectrum. However, it is very difficult to 
establish Hagedorn's conjecture in a compelling way from the empirical data. In part this 
is a practical issue; one would need to extract hadron masses for hadrons up to much larger 
masses to get compelling evidence for an exponential growth. Moreover, underlying this 
practical issue is an important theoretical one: highly excited hadrons are not particles; 
they are resonances and as such do not have well-defined masses. The mass parameters 
can only be extracted from partial wave analysis of various scattering processes using some 
model dependent assumptions. Such model-dependence is quite weak for well-isolated nar- 
row resonances and for these one can state masses with some level of confidence. However 
as resonances in some channel become wide or close to each other, such model dependence 
grows and it becomes difficult to isolate resonant state in a meaningful way. Moreover, 
any model dependence in the meaning of a hadron's mass makes the issue of the density 
of hadronic states intellectually problematic. 

Before attempting to deal with the problem of ill-defined hadron masses, it is useful 
to understand why one might expect QCD to have a Hagedorn spectrum. Recall that 
Hagedorn spectra arises automatically in simple string theories [7] with unbreakable and 
noninteracting strings. It is noteworthy that string theory was originally formulated as 
a theory of strong interaction. Moreover, given confinement it is plausible that highly 
excited states in QCD should act stringy. For the case of pure gauge theory there is strong 
evidence [8] that widely separated static quark sources have a linearly rising potential, (i.e., 
an area law for the Wilson loop). This arises because for widely separated sources, the 
flux arranges itself into tubes with a characteristic width and fixed energy per unit length 
[8]. It is plausible that for highly excited states, which would be expected to have flux 
tubes which are much longer than their width would act dynamically as strings and as an 
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effective string theory would naturally give rise to a Hagedorn spectrum. Mesons in such a 
picture are interpreted as open strings. However, this picture is flawed. It is based on pure 
gauge theory, in which confinement implies unbreakable strings. In QCD, with dynamical 
quarks, flux tubes can break. Indeed, this is the same issue as noted above — the fact that 
flux tubes break implies that mesons decay and thus can only be seen as resonances with 
nonzero widths. 

It is not clear whether there is a clean way to deal with this issue in an unambiguous 
way for QCD in the physical world. However, if one focuses on the large N c limit of QCD 
[9, 10], the issue vanishes. As the large N c limit is approached, meson decays are suppressed 
by a factor of 1/N C ; flux tubes do not break and mesons become stable. The goal of the 
present work is to show that at large N c QCD must have a phase transition. This is at 
least a well-posed theoretical question. Of course, the question of whether or not this is 
of phenomenological relevance depends on how close the N c = 3 world is to the large iV c 
world. 

The first derivation of a Hagedorn spectrum in some variant of large N c QCD was 
done by Kogan and Zhitnitsky [11], who explicitly computed the spectrum of large N c 
QCD in 1+1 dimensions with adjoint fermions and showed that it possess a Hagedorn-type 
behavior. Ideally one could similarly compute the spectrum for large N c QCD for more 
than one spatial dimension. However, in practice we do not know how to solve for the 
spectrum. Numerical studies using lattice QCD are poorly suited for extracting high-lying 
stars. There was a study of the large N c glueball spectrum based on a numerical treatment 
of a transverse lattice QCD in a light cone formalism [12]; while the results are consistent 
with the Hagedorn spectrum, the evidence was not definitive. There are, however, indirect 
ways to probe the issue. One is by the study of QCD thermodynamics. It is well known that 
large N c QCD has a first order phase transition to a quark-gluon plasma phase [13, 14] with 
the latent heat growing as N% . This transition tells us nothing about a Hagedorn spectrum. 
However, systems with first order transitions can superheat and thus a hadronic phase can 
exist about T c . The Hagedron spectrum and the noninteracting nature of hadrons at large 
N c implies a maximum temperature for this superheated hadronic phase [4, 15]. Moreover, 
it is practical to study this superheated phase in lattice QCD for moderately large N c and 
thus get indirect evidence for a Hagedorn spectrum. 

An alternative indirect way to demonstrate a Hagedorn spectrum for large N c QCD for 
3+1 dimensions was outlined in ref. [16]. The argument relies only standard and generally 
accepted properties of QCD. Confinement in its basic sense that all physical states are 
color singlets plays a critical role as does asymptotic freedom. In addition, the approach 
requires some plausible assumptions about the validity of perturbation theory to describe 
the correlation functions at short times. However, the approach explicitly assumes neither 
that the hadron dynamics is stringy in nature nor that the confinement is manifest through 
an unbroken center symmetry. The argument relies critically on the fact that the number 
of independent local operators with given set of quantum numbers grows exponentially 
with the mass dimension of operators. This approach is similar in spirit to the ideas of 
Kogan and Zhitnitski [11]; it also has elements which are reminiscent of refs. [17-19]. 

The principal purpose of this paper is to generalize the argument of ref. [16]. The 
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version of the argument in ref. [16] does not apply to 2+1 dimensions; here we develop a 
variant of the argument which is applicable to both 3+1 and 2+1 cases. We also simplify 
and clarify the arguments in ref. [16] and improve it in significant ways. The principal 
improvement is in the treatment of perturbative corrections to correlations in section 7 of 
this work. In ref. [16] corrections due to certain classes of diagram were shown to have 
a required behavior and it was suggested that the general case ought to behave similarly. 
Here a complete demonstration that this is true is given. 

2 Outline of procedure 

We start by noting that at large N c meson widths go to zero. Thus the spectrum of mesons 
is unambiguously defined. To begin the analysis we define two functions characterizing the 
spectrum of hadrons. N(m) is defined as the number of hadrons with mass less than m, 
and W(m) defined as the sum of the masses of respective particles 

W(m) = ^mi = mN(m) - j dfiN(p) . (2.1) 

o 

It is easy to see that if one of them grows exponentially so does the other. 

Next, we explicitly construct a sequence of sets of local operators with fixed mass 
dimension (labeled with n) which grows exponentially. The number of operators in each 
set of the sequence is 

N = A n . (2.2) 
Our goal will be to demonstrate that at sufficiently large n the following inequalities hold 

N(an + b) > V ; (2.3) 
V > W{m N ) , (2.4) 

where V is the negative logarithmic derivative of the trace of a certain matrix of correlators 
and a and b are constants with dimensions of mass. The key feature is that the left-hand 
side grows linearly with n and no faster. 

If we now assume that the number of hadrons is bounded from above N(m) < exp(am) 
we can easily derive from eqs. (2.1)- (2.4) following expression: 

aalog A (e)m + b > m . (2-5) 

a 

As m — > oo there is a contradiction unless a > , 1 , s and the assumption that function 

— alog A (e) ^ 

N(m) is bound by exponential is false. Consistency requires 

N(m) > exp ( - — - m) . (2.6) 
Valog A (e) J 

Consequently we obtain the Hagedorn spectrum. 

Inequalities (2.3) and (2.4) are somewhat subtle and their derivation will be described 
in detail in the following sections. 
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3 Sets of local operators 

To proceed further we need to construct sets of composite single-color-trace color-singlet 
local operators. The matrix of correlators of operators in these sets is the core of the 
argument. The single color trace nature of these operators will ensure that at large N c 
each of these operators when acting on the vacuum makes a single hadron (provided that 
one assumes confinement) [9]. These operators need to have the property that at large 
N c , the correlator between two distinct operators in the set must vanish as the distance 
between the operators goes to zero. For simplicity we consider operators which transform 
as Lorentz scalars; these are guaranteed to produce spinless hadrons when acting on the 
vacuum (for any spatial dimension) . It is sufficient to show that that the number of spinless 
hadrons grows exponentially to establish a Hagedorn spectrum. 

The operators we use in our construction need to be different for 2+1 and 3+1 dimen- 
sional QCD. We will construct our operators out of some basic building blocks. In 3+1 
dimensions, these building blocks are the following two types of operators: 

01 = const • F^F^ , 2 = const • F^F^ , (3.1) 

and in 2+1 dimensions: 

Oi = const • F a/3 F^ F al p,F al P' F a „ p „F a "P" , 

(3.2) 

where the constants may be chosen for convenience and do not affect any results. These 
operators are not traced over color. Thus, in the large iV c limit these become pure color 
adjoint operators. It is easy to see that such operators are linearly independent — one is a 
scalar and one is a pseudoscalar. 

From these basic building blocks, we create the individual color-singlet operators in 
the following way: 

Jh,l 2 ,...,l n = qO h Oi 2 . . . O in q , (3.3) 

where Is are either 1 or 2 and n is the total number of the operators O inserted between 
quarks. We construct sets of these operators each of which has the same value of n. Thus, 

Si 

So 

s 3 

Ultimately we consider these a sequence of sets where the n th element of the sequence is 
S n . We will then focus on the behavior as n becomes large. By construction, the number 
of currents in the n th step of the sequence equals N = 2 n . Additionally, all currents in a 
given set S n have the same (naive) mass dimension, 4n + 3 for 3+1 dimensions and 12n + 3 
for 2+1 dimensions, respectively. To sum up, we created a sequence of sets of operators 
where the number of elements in each step grows exponentially with a mass dimension. 



const 



7/i v 



F^'u" 



= {Ji,J2\ = {qOiq,q0 2 q} , 

= {Jll, J12, -hi, J22} = {qOiOiq^O^q, ■ ■ ■ } , 

= {yiii, ^112, J\2U J122, J211, J212, ■ ■ ■ } , 
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4 Matrix of current correlators 



Let us define a sequence of correlator matrices Il( n ) between two space-time points. Their 
matrix elements read 

n^(x-y) = (ji(x)J b (y)) , (4.1) 

where currents J a ^ G S n . The dimension of such matrix is equal to the number of currents 
in the respective set, i.e., 2 n . To keep things clear, we will be using the following notation: 
matrices elements will always be written with explicit indices whereas matrices themselves 
will be indicated by boldface (n a f, II). 

The large N c limit plus the assumption of confinement guarantees that every current 
generates only single meson states; the widths go to zero as N c approaches infinity. Thus, 
the spectral decomposition of the correlator is given by 

Ut,W) =E/ IM ■ ("> 

Here, c a k is the amplitude that the current a creates the particle k. 

Using eq. (4.2) we can write the matrix II (without loss of generality we can take 
y = 0) 

n^M) = E / (§3 c <** > ( 4 - 3 ) 

where = c* ak Cbk is the matrix of coefficients and A(i,x;mfc) is the propagator for 

a noninteracting scalar of mass m\- . Our matrix can be viewed as the Kallen-Lehmann 
spectral representation with the spectral function p proportional to Dirac delta functions. 
It is a straightforward consequence of the large N c [9, 10] limit, planarity of diagrams, and 
confinement. 

We will study only a correlation in time (x = 0) and perform an analytic continuation 
to an imaginary time (r = it), thus 

n^V) = E / (0 c *j> A (w) • ( 4 - 4 ) 

5 The relation between matrix of correlators and masses of mesons 

In this section we derive the eq. (2.4) which is at the heart of the demonstration of a 
Hagedorn spectrum. We note that a derivation of this was given in ref. [16] . Unfortunately 
that derivation contained an error. The result, however, is correct and a valid derivation 
is given here. 

Before proceeding, it is useful to recall that it is standard to use current-current cor- 
relation functions in lattice QCD to extract the lattice hadronic state in a given channel 
[8]. The following relations are used in this context: 

l im _Aio g (J( r )J(0)) = m , (5.1) 

t— s-oo or 

-Alog<J(r)J(0)> > m , (5.2) 
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where mo is the lowest mass state. The goal here is to generalize the relation in eq. (5.2) 
to the case of a matrix of correlators Il( n ) in large N c QCD. 
Prom eq. (4.4), the matrix elements are 

ni^r) = (4(r)J 6 (0)) = E / (§3 C a b,k ^ , (5.3) 

where a, b = 1 ... 2 n . We study the following expression and will prove that it is greater 
than the sum of the lowest 2 n masses with scalar quantum numbers. 

= -^-TrlogllW = Tr J-ri^nW" 1 ) , (5.4) 

where the dot stands for the derivative with respect to Euclidean time r. 
First, let us define the following quantity: 

compared to the derivative with respect to r: 

" ni? (r) = £ / U C ab ,k ±- > (5-6) 

k 

One can easily prove that the trace (5.4) with n instead of II is always smaller than the 
trace with the dotted matrix: 

Tr (-n) n _1 > Tr (-11) II 1 , 

E(^c| - nmK 1 > X>d - n|^ c )A- 1 , (5.7) 

c c 

where \ip c ) are eigenvectors of matrix II, and A c are corresponding eigenvalues. All eigen- 
values, A c , are positive since II is a positive definite matrix. The left part of inequality 
(5.7) represents the average of energy (which includes momenta) in certain states, which 
is always greater than an average of the analogous masses, that correspond to the term on 
the right-hand side. The average is taken with positive weights A" 1 so the inequality holds 
overall. 

Next, we split the matrix II into two parts, A and B, where A couples only to the 
first 2 n masses and B to the rest. 

d 3 P ^ 1 -E k r 

e 



A ab i T ) = JL J (2^)3 Cab > k ^ 



00 f' r\3 ~" 1 



k=2 n +l 



One can easily see that such splitting can be done for the matrix itself, II = A + B, its 
derivative, n = A + B, as well as for the matrix, n = A + B. 
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Our goal is to show that we can relate the II matrix to the meson masses. We start 
with the following simple matrix identity (which is valid for every n): 



Tr ((-A - B) (A + B)" 1 ) = Tr ( - 1 VX ( -L(-A)-^) VX , 1 
V V M / VVA + B Vv^ WaJ Va + b 



+ Tr ( -=L= VB f -±=("6)4=1 V5- 



VA + B Vv 7 ^ WBJ VA + b 



(5.9) 



Next we introduce the definitions 



A=-Jr(-A)-Jr, B=i(-fi)>. (5.10) 



which allows us to recast eq. (5.9) in the following way: 

Tr C(-A-B) (A + B)" 1 ) = Tr (x f AX + Y f B y) . (5.11) 



The relation between X and Y operators X^X + Y^Y = 1 guarantees that they can be 
simultaneously diagonalized and their eigenvalues are related; 

X+Xhfc) = Pi m , YtY|^i) = (l-Pi)|Vi), (5-12) 

with < pi < 1. We can see that the states defined as |x?) = l/^Xl^j) forms an 
orthonormal basis of states, as do the states |Tj) = — p% Y\ipi). Thus, we can rewrite 

the trace into the following form. 

Tr (x+AX + Y+BY) = ^(^IX+AX + Y+BY|^i) 

i 

= ^(Xi|A|xi> + (l-Pi)«Ti|B|T i )-(Xi|A|xi)) (5-13) 

i 

> TrA = Tr— - . (5.14) 
A 

The second term of the right-hand side of eq. (5.13) is always positive since the biggest 
possible eigenvalue of operator A is less than or equal to rri2 n (where rrij is the mass of the 
j th state) , while the smallest eigenvalue of B is always greater or equal to m^+i- This is 
a straightforward consequence of eq. (5.5), (5.8), and (5.10). 

The currents entering operator A can be reorganized in such a way that the first 
current couples only to the first meson, second current only to the second meson, etc. It is 
obvious that this process will lead to the diagonal matrix whose every element is equal to 
the mass of the n th particle. Consequently 

-AW ^ 



- 8 - 



The quantity on the right-hand side is exactly the function W(ni2^) defined by eq. (2.1). 
Consequently, we derived the inequality of (2.4) required for the establishing of a Hagedorn 
spectrum. 

on 

V = -— Trlogn> y^m k = W(m 2 n) . (5.16) 
dr 

k 

6 Matrix of correlators in an asymptotically free regime 

In this section, we derive the condition (2.3) stating that the trace of the logarithm of a 
corellator matrix grows at most linear in n (where n labels the step in the sequence). 

For sufficiently small times, asymptotic freedom allows us to treat fields inside the 
currents as non-interacting, so we can decompose the correlator to a product of single- 
particle propagator functions. Additionally, the large iV c limit, in which we are working, 
guarantees that the whole matrix of correlators is diagonal. Doing the trace, we obtain 
2 n terms with the same structure. Consequently, we can focus solely on the one current- 
current correlator and its logarithmic derivative, and investigate how it grows with n. 

In the asymptotically free region, the structure of the overall correlator is given simply 
by the dimensional analysis. Since the time scale r is the only dimensional parameter left, 
and our current J has the mass dimension 4n + 3 for 3+1 dimensions, and 12n + 3 for 2+1 
dimensions, the right-hand side equals 8n + 6, and 24n + 6, respectively. Thus 



jj(n) I <5 a b const t 8n+6 for 3 + 1 dimensions 
ab &ab const r~ 24n+6 for 2 + 1 dimensions 



and the trace of logarithmic derivative equals 

-A T rlo £ nH-/( 2n )^ for3+ldimensions 

^ log - | (2n) 24^6 for 2 + 1 dimensions 

If we assume that the matrix of correlators is effectively at its asymptotic value up to some 
small corrections for r < tq with tq independent of n, we reproduced exactly the inequality 
condition (2.3) necessary for establishing the Hagedorn spectrum. 

From the argument in section 2, the preceding implies a Hagedorn spectrum where the 
value of Hagedorn temperature corresponding to our sets of currents is 

f 81og 2 (e) 



H j 24iog 2 (e) £ Qr 2 _|_ i dimensions ^ ^ 

At this stage, we have shown that a Hagedorn spectrum emerges in the QCD in 3+1 
and 2+1 dimensions if a certain assumption is met — namely that correlators are to good 
approximation at their asymptotically free value for r < tq for all n. 

7 Perturbative corrections 

In the previous section, we neglected all possible interactions between gluons. Such an 
assumption was justified by the fact that the QCD is in asymptotically free regime for 
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short times. Provided that this condition is met, we have proved QCD has a Hagedorn 
spectrum. The critical question is then the circumstances for which there are no large 
corrections to the asymptotically free result. 

The standard way to include the effects of interactions for the correlators at short 
times is via perturbation theory. However, the region where perturbation theory is valid is 
certainly limited; as the time increases perturbative corrections grow and ultimately push 
the system outside the region of validity of perturbation theory. Here we will rely on the 
standard assumption that perturbation theory accurately describes correlation functions 
provided that they are small. That is, in the region where pertrubative corrections are small 
they will dominate over all nonperturbative effects. We note that this is not a rigorous 
mathematical theorem but it is the basis of standard analysis of QCD correlation functions. 

Given this assumption, the critical issue we need to address is how do perturbative 
corrections scale with n? The goal is to show that at fixed r perturbative corrections, at any 
fixed order, to the quantity 2~ n ^:Tr logn^ scales at most linearly with n for any fixed r 
(in order to satisfy the inequality (2.3)). If we can demonstrate this, we have demonstrated 
a Hagedorn spectrum given the assumptions stated above. To see why, imagine doing the 
following calculation: start at some fixed but large n and compute 2 -n ^Tr log n^") in 
perturbation theory to some order. Next decrease the value of r so that perturbative 
corrections are sufficiently small that the quantity is close to its asymptotically free value, 
up to corrections which are a small fraction of the total. It is always possible to find a 
value of r for which this is true since the system becomes asymptotically free as r — > 0. By 
assumption this is the regime in which perturbation theory is trustworthy. Having done 
this at some fixed value of n, one next increases n keeping r fixed. Since the perturbative 
corrections to 2 -n g^Trlogn ( ' n ) have been demonstrated to scale at most as n and since, 
as seen in eq. (6.2), the leading behavior from the asymptotically free region also scales 
with n, the fraction size of the correction is independent of n. Since the corrections to 
the asymptotically free result were small at the original n, they remain small at all n 
including in the limit of n — > oo. This is sufficient to show a Hagedorn spectrum, given the 
assumptions stated above, and given the result of section 6. 

In the remainder of this section we show that perturbative corrections to 
2~ n ^:Tr log n( n ) do, in fact, scale with n no faster than linearly and thereby complete 
the demonstration. We note that an argument that this is the case was presented in 
ref. [16]. In that work, the results of certain classes of diagrams were presented and shown 
to be consistent with the needed result. It is was suggested that the structure of the 
quantity ought to ensure that result continued to hold for all classes of diagram, however, 
no demonstration of that was given. Here we construct a general argument why the result 
will hold for all classes of diagram. 

In order to get more insight into the correlation functions, we first look at the contri- 
bution from non-interacting gluons such as in the Feynman diagram illustrated in figure 2. 
It will allow us to extract some properties that we will generalize later when we include 
interactions. In the non-interacting case, the correlator matrix is diagonal due to large iV c 
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n 



Figure 2. The Fcynman diagram for a correlator in the asymptotically free regime. The blobs 
indicate the currents. 

limit. Each propagating gluonic operator 0\ inside the current contributes the same way 

%i,Jj) = (Oi{T) O,(0)) free = Sij 4 ee (r) , (7.1) 

where the superscript 1 indicates that this represents the propagation of a single gluonic 
operator. The 5ij is a result of the choice of operators (3.1), (3.2) which were picked 
precisely because of this property. The fact that the free correlator is the same for both 
operators is a result of dimensional analysis which requires them to be proportional to each 
other; a choice of the constants defining the operators can fix the proportionality constant 
to unity. We want to emphasize that this object is not gauge invariant, however, it will 
be only a part of the overall correlator which will be gauge invariant. The small letter tt 
indicates that the quantity is already a function, not a matrix. Such convention will be 
applied also in the following text. Let us also remind the reader that the matrix elements 
explicitly indicated indices, whereas the analogous matrices are denoted in boldface. 

The correlation function of the whole current consists of n internal lines and is bounded 
by two quark lines. Thus its matrix elements are given by 

ni? (r) = 6 ab nlJr) [^(r)] B 4 ee (r) . (7.2) 

where ^ tcb (t) is the free quark propagator traced over Dirac indices. Recall that the 
indexes a, b goes from 1 to 2" . The quantity of interest is the derivative of trace of the 
logarithm, 

2-" A T r lognW = 2^1og7r| ee + nAl og7r i ree . (7.3) 

The first term is independent of n (and comes for 2 quarks on the boundaries), while 
the second one is linearly proportional to n. This result is in agreement with the results 
obtained in the previous section. 

Diagrams with interactions are more complicated. Let us work in a case, when all 
interaction up to order a 1 are included. First, we restrict our attention to the case where 
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n 



Figure 3. An example of a Feynman diagram where interactions do not couple distinct gluon lines 
connected to the sources. 

all of the interactions act on single gluons that are connected to the sources (see figure 3). 
While more than one of these gluons may be involved, there are no interactions which 
couple distinct gluons coupled to the source in this class. Effectively, these interactions 
lead to the modification of a free propagator similar to the contribution of self-energy 
correction. 

4ee(0 -> TTfreeW (1 + c(t)) , (7.4) 

where (1 + c(t)) is a pertrubative correction. Note that c depends on the order to which 
we work in perturbation theory, but it is well defined at any given order. Such correction 
can appear on any internal gluonic line, so that if these were the only types of diagrams 
contributing we could write the total correlator as 

nS?(r) = 5 ab 7r? ree (T) [4cc(t) (1 + c(r))] n vr| ee (r) . (7.5) 

This structure actually contains more information (higher order in a) than necessary, but 
certainly contains all combinations that are required to order a 1 . The trace of the logarithm 
now reads 

2~ n log nW = 2 i- log 7r<? ee + n A i g ^ + n _l i og (i + c ) . ( 7 . 6 ) 

Although we have yet to calculate c, it is clearly independent of n. So the total expression 
grows, again, at most linearly with n. Of course, this result is wrong — the class of diagrams 
we considered was chosen artificially and does not correspond to the full perturbative result 
at any order in a. However, it does illustrate how the logarithmic structure combined 
with factorizing point-to-point correlators for the individual lines yields total perturbative 
corrections which grow at most with n. Our main task is to show that a structure enriched 
with inter-gluonic interactions obeys the same at most linear n dependant growth. 

Up to now, we have ignored possible interactions between gluons. In order to deal with 
them in a simple way, we artificially divide the total number of internal gluon lines n into 
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n1 



n1 



n2 



Figure 4. Division of gluon lines into 712 clusters each containing n\ lines. 



n,2 clusters each containing n\ gluon lines, that is n = n\U2- We will impose the conditions 
n\ S> 1, fi2^> 1; and rii ^> the order in perturbation theory to which we are working. For 
the beginning, let us assume that all interactions occur within clusters, as is illustrated in 
the figure 4, i.e., there is no internal line between clusters or between these clusters with 
quark lines bounding them. Obviously, in doing this we neglect a certain class of diagrams 
for now. However, we will subsequently show that correction due to their inclusion does not 
affect the leading n behavior. The reason for this is that for a fixed order of perturbation 
theory only a very small fraction of all Feynman diagrams will connect different clusters 
and this fraction is small enough to alter leading behavior. 
The propagation within one cluster is given by 



where [1 + C'(r, ni)] represents the effect of interactions within one cluster. Note that 
II C is a matrix of the dimension 2 ni x 2 m . The C'(r, n\) depends on the order to which 
we work in perturbation theory, and, obviously, it depends on the size of the cluster, the 
number of internal lines n\. 

Now, we need to define a mapping from the "cluster" space, which has the dimension 
2™ 1 matrix to the "overall correlator" space with the dimension 2 n , D(M). More precisely, 
we actually need ni different mappings D^ fc )(M). each corresponding to a different cluster. 

(k) 

Our mapping will be a 2" x 2 n matrix D b . Note that indexes a and b can be repre- 
sented by a sequence of n numbers, a = (aid? ■ ■ ■ a n ), b = (6162 • • • b n ), where cjj, 6j = 0, 1. 
It is a straightforward consequence of our original definition of currents (3.3), each of them 
being constructed from n building block operators (translated into n internal gluonic lines 
in each diagram). Recall that we divided n lines into r&2 clusters of n\ elements . Moreover 
since we neglect interactions between clusters, we want to treat them independently. Thus 
it is useful to define a mapping D( fc ) to work in such a way that the first one, DW, affects 
only first n\ subindices within indices a and b, the mapping corresponding to the second 
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one, D^ 2 ), affects pieces n\ + 1 to 2ni, etc. Analogously, all n 2 mappings corresponding to 
all possible n 2 clusters are denned. Specifically, 



D ^ 1 a 2 ...a n )(b 1 b 2 ...b n )( M ) ~ M (ai...o„ 1 )(&i...&„ 1 ) ]J 6a i b ' ' 

l= ni +l 

ni n 

D (a 1 a 2 ...a n )(b 1 b 2 ...b n )( M ^ = M (a ni+1 ...a 2ni )(b ni + 1 ...b 2ni ) Yl^<H b t II ^ a l' b l' ' 

1=1 l'=2n 1 +l 

and the general form 

(fc-l)ni n 

D( ai a 2 ...a n )(b 1 b 2 ...b n )(^ / ty = ^(a{k-i)n 1 +i---akn 1 )(b(h-i)n 1 +i---bk-n. 1 ) II <Wl II <W&j/ ' 

Z=l Z'=fcni+1 

(7.8) 

where M is an ni x ni matrix. Note that such matrices are a straightforward generalization 
of the previous case where the equivalent of Il c was just number (matrix lxl) and the 
matrix D was diagonal. It is worth mentioning that the 2 n dimensional overall space 
can be factorized as a product of n 2 subspaces with dimensions 2™ 1 with each subspace 
corresponding to one particular cluster. So, the vectors have the form 

\V) = I^W) <g> |?/ 2 )) <g> • • • ® |^ n2 )) , (7.9) 

and the spirit of the mapping D( fc ) is 

D^(M) = 1« <g> 1< 2 > ® • • • ® ® ■ ■ ■ ® 1^ . (7.10) 

From the construction of D^^(M) it is easy to show that Tr logD^ fe ^(M) = Tr logM. 

Using this notation and imposing the condition that we neglect all diagrams connecting 
clusters, the overall correlator contains a product of n 2 matrices corresponding to the 
respective clusters 

n(n) = <ce( T ) D«(n c ) x D( 2 )(n c ) ... n^\n c ) ^(r) . (7.11) 

Using the general property that Tr log(AB) = Tr log A + Tr logB and the property that 
Tr logD^^M) = Tr logM, it is straightfoward to show that 

2~ n A Tr lognW = 2 A logvr| cc + m« 2 A log4 ee + n 2 A log(l + c'(m)) (7.12) 

where c'(ni) = exp (Tr log (1 + C')) — 1. 

The first two terms were already discussed after eq. (7.3). The second term becomes 
obvious once one recalls that mn 2 = n. The third term requires more care. We can denote 
the expression ^log(l + c'(ni)) as f{n\) , i.e., some at present unknown function of n\. 
However, our choice of clusters is completely arbitrary. Provided our assumption that 
diagrams connecting clusters does not affect the leading behavior is correct, we can switch 
n\ and n 2 and the result must remain unchanged. This requires nif{n\) = ^1/(^2) an d 
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n1 n1 
n2 

Figure 5. An example of a Feynman diagram with interaction between two clusters. 

consequently f(m) must be a linear function of m. Thus, the third term on the right-hand 
side of eq. (7.12) is also proportional to n = n\n%, just as in eq. (7.6). 

To complete the demonstration we need to show that the inclusion of diagrams with 
interactions connecting individual clusters does not affect the leading scaling. An example 
of such diagram is in figure 5. These effects can be accounted for by including a matrix 
(1 + C") between the matrices corresponding to clusters. Additionally, one should include 
the interaction between the quarks on the boundary and the neighboring clusters of gluons 
(1 + C q ). The total correlator matrix reads 

nW(r) = < ree (r) (1 + C 9 ) x D^(n c ) x (1 + C") x D^(n c ) ... 

(1 + C") x T>^(U C ) x D«(II C ) (1 + C 9 ) ttIJt) . (7.13) 

The key point here is that the matrices C" and C q must be independent of n\. The reason 
for this is that by construction, n\ is much larger than the order in perturbation theory 
to which we are working. At large N c only planar diagrams contribute. Thus if we are 
working at order a! s , a diagram connecting two clusters can at most go I — 1 gluon lines into 
the cluster. Since this is less than rt\ it cannot go across the cluster. Thus, the dynamics 
in C" does not know how large the cluster is and must be independent of n\. It is clearly 
independent of rii either. 

Using an analogous argument as earlier, the trace of the logarithm reads 

2^ A Tr log OW = 2 A log (vr | ce(1 + ) + nin2 ±. log ^ 

+n 2 A log(l + c'(m)) + (n 2 - 1) A log(l + c") (7.14) 

where c q and c" are defined analogously to d . The first term does not scale with n; 
the second term is directly proportional to n. Defining f(n\) = log(l + d{n\)) and 
g = £p log(l + c"), the last two terms can be written as n2(f(n{) + g) — g. Again exploiting 
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the fact that we can switch m and n 2 without affecting the result (provided the order in 
perturbation theory is less than both ri\ and n 2 ) we obtain the consistency condition that 

n 2 {f{n l )+g)=n l {f{n 2 )+g) (7.15) 

which yields f(m) = bm — g where b is a constant. Thus, the effect of including the 
interactions between the clusters simply fixes the subleading behavior in /. Taking the 
last two terms together and exploiting the fact that n = n\n 2 we have bn — g which grows 
linearly in n. Consequently, the right-hand side of eq. (7.14) grows at most linearly with n, 
and the inequality (2.3) is satisfied if the interactions are included via perturbation theory 
at any fixed order. 

With this we have completed our demonstration that a Hagedorn spectrum arises in 
large N c QCD in both 2+1 and 3+1 dimensions. This demonstration depends on one crit- 
ical assumption: that perturbation theory accurately describes the trace of the logarithm 
of a matrix of point-to-point correlation functions in the regime where the perturbative 
corrections to the asymptotically free value are small. 
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